Novel classical ground state of a many body system in arbitrary dimensions 

G. Date*, Pijush K. Ghosh*, and M. V. N. Murthy* 
The Institute of Mathematical Sciences, Madras 600 113, India. 
(February 1, 2008) 



OO 
r—i 

CU- 
in 



> 

O 

m 

(N 
O 
OO 

o\ 

1 



O 
o 



X 



The classical ground state of a D- dimensional many body 
system with two and three body interactions is studied as 
a function of the strength of the three body interaction. We 
prove exactly that beyond a critical strength of the three body 
interaction, the classical ground state of the system is one in 
which all the particles are on a line. The positions of the 
particles in this string configuration are uniquely determined 
by the zeros of the Hermite polynomials. 

PACS numbers: 03,20.+i, 03.65. Ge, 24.10.Cn 



In this letter we consider the following Hamiltonian 
which displays remarkable classical properties. The clas- 
sical Hamiltonian of the model for N interacting parti- 
cles, in D dimensions, is, 
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where gi and 172 are in general arbitrary but positive cou- 
pling strengths of the two and three-body interactions. 
We have set the mass m = 1, oscillator frequency u> = 1 
and fij =fi — fj. 

We prove exactly that when 52 > <?i the classical 
ground state of this Hamiltonian is unique (modulo ro- 
tations) and is described by all the particles along a line 
with locations determined by the zeros of the Hermite 
polynomial H^. 

There are several reasons to consider this Hamiltonian. 
The model Hamiltonian given above has a very interest- 
ing limit. When 



gi =g(g + D-2) 
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it has been shown that the quantum mechanical ground 
state and an infinite tower of excited states may be ob- 
tained exactly in D-dimensions jjJD. ( I n D=2, some of 
these solutions may be obtained from the Hamiltonian 
of the anyons in two dimensions in the zero total angu- 
lar momentum sector In particular for the ground 
state solving the Schroedinger equation HijiQ = Eq^/jq we 
obtain, 
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and the quantum ground state energy is given by, 



E = -[ND+gN(N-l)], 
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where g > 0. In the limit g — ► 0, the system describes 
N free bosons confined in a harmonic oscillator potential 
without the centrifugal barrier. 

For D = 1, the model is the celebrated Calogero- 
Sutherland model (CSM) 0. The three-body term is 
identically zero in this case. It is known that the classi- 
cal ground state of CSM is the one in which the location 
of the ./V particles are determined by the zeros of the Her- 
mite polynomial Hn ||- The CSM is of course exactly 
solvable quantum mechanically. 

For D = 2, when 172 = 0, the model reduces to a model 
with only two-body inverse square interaction which has 
some interesting applications to quantum dots and is con- 
sidered in detail in ref. ||. Unlike in one dimension, 
however, the three body term does not vanish in higher 
dimensions and plays a crucial role in the analysis that 
is to follow. Similar to the one dimensional CSM the 
norm of the ground state wave function in eq.(^), when 
g 2 = gi = g 2 , can be identified with the joint probabil- 
ity density function of the eigenvalues of complex ran- 
dom matrices for a particular value of g Also, the 
ground state wave function is exactly found to be the 
gauge transformed Laughlin's wave function and exhibits 
off-diagonal long range order ^|. For arbitrary g\,g 2 , 
however, the quantum problem may not be solvable even 
for the ground state. 

We also note that the three-body interactions have 
been used in several areas like nuclear, atomic and con- 
densed matter physics. The main reason usually is to 
account for the many-body effects in systems of particles 
or clusters. It is also known that the three-body interac- 
tions may be used to bring about structural transitions 
in the classical ground state. Recently, Weber and Still- 
inger |J] have proposed a form of the three-body inter- 
action for inducing structural transitions in the classical 
ground state of many particle systems in their molecular 
dynamics simulation. The three-body interaction used 
here bears a resemblance to the one proposed by Weber 
and Stillinger. 

In what follows we keep g\ and gi arbitrary but positive 
and analyse the classical ground state of the Hamiltonian 
in eq.(|lJ) in D-dimensions. Defining 
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the Hamiltonian can be written in the form, 
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Therefore A = corresponds to the special case when 
gi = <?2- Here g may be chosen to be either positive or 
negative. The results that follow, however, cannot de- 
pendent on this choice since the Hamiltonian and the 
equilibrium conditions depend only on g 2 . We choose to 
work with g > since the analysis is simpler. The classi- 
cal equilibrium configurations (fixed points) are obtained 
by solving, 



Pk = 0, 
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where k — l,...,N and the derivatives are taken with 
respect to configuration coordinates ffe. The energy at 
any of these extrema are given by, 
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where fi are solutions of the eq. ( jlpj ) . The last part of the 
above equation follows from the fact that the solutions 
of the equilibrium equations imply, in general, 
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Thus the interaction energy is equal to the confinement 
energy for all equilibrium configurations. 

Before analysing the Hamiltonian further we first es- 
tablish properties of the three-body potential V3. Sym- 
metrising V3 explicitly, we have 
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and using the fact that + r\/& + r*fci = 0, it is easy to 
show that 
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where ^-fc is the angle between fj and r}/.. From eq. (J14|) 
following results follow independent of whether f are so- 
lutions of the equilibrium equations or not: 

1. Clearly V 3 > for all f. 

2. Further V3 = means that each term in the sum- 
mation must be zero since each term is separately 
a square. Therefore, dijk must be either zero or n. 
This means that all particles must be on line when 
V3 attains its global minimum. 

3. By a similar calculation it is easy to see that, 
Vil/3 = for all i if all the particles are on a line, 
that is, f = Sii, where the unit vector e denotes 
an arbitrary direction in the D dimensional config- 
uration space. 

4. Therefore, if f = Sje is an extremum then the solu- 
tions of eq. (0) as well as the energy in eq. ([ll]) are 
independent of A. Furthermore this energy attains 
its absolute minimum value namely, 
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To compare qualitative features for A positive, nega- 
tive or zero we consider special configurations such that 
all the particles lie on a plane passing through the ori- 
gin. The plane in fact can be chosen to be the x-y plane 
(say) because of the overall rotational invariance. The 
equilibrium configurations among these are governed by 
the same equations. Of these planer configurations we 
consider two special ones which are also solutions of the 
equilibrium equations, namely, (1) all N particles on a 
circle with centre at the origin and (2) (N-l) particles on 
a circle and one at the centre. For these configurations 
the respective energies are given by, 



E, 



O 



2A N -2 
Sg 2 N -V 



Eq=E, 



line 




2A (N + i)(N-3) 



N 2 



(16) 



(17) 



where Eu ne is given in eq.(|15|) and refers to the energy 
of the line configuration (which is independent of A). 

(a) A < 0: For this case the line configuration can not 
attain the absolute minimum of the energy since the spe- 
cial configurations mentioned above have lower energy. 
For arbitrary N the ground state is indeed hard to find 
though it is clear from the form of V3 that these are nec- 
essarily non one dimensional configurations. Among the 
sub-class of planer configurations, or equivalently D = 2, 
we know that for N < 5 circle configuration has the low- 
est energy where as for 6 < N < 8 circle-dot has the 
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lowest energy. For N > 9 multishell configurations have 
the lowest energy Jk|. For D > 2, of course these in 
general are not classical ground states. Note that even 
though A < 0, by definition A > — g 2 since the strength of 
the three body term is necessarily positive. This ensures 
that the energy is real in these configurations. 

(b) A = 0: The absolute minimum of energy is attained 
for all configurations satisfying the equation Ri — 0. The 
energy in all these configurations is given by, 



E = g 
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In this case the line, the circle and the circle-dot are 
all degenerate. Within the sub-class of planer configura- 
tion (equivalently for D — 2) numerical simulations show 
that there are many other configurations which are also 
degenerate. For D = 2, the quantum mechanical ground 
state can be obtained exactly in this limit as noted at 
the beginning. The classical energy is the same as the 
quantum mechanical ground state energy apart from the 
contribution from the zero point fluctuation. Note that 
this is rather special for D = 2. 

(c)A > 0: The line configuration is certainly a ground 
state since the absolute minimum of the energy is at- 
tained. Because V3 is positive definite or zero, any other 
extremum must necessarily have higher energy, 
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Thus for all A > 0, the classical ground states are line 
configurations with f*j = s^e, where e defines an arbitrary 
direction in the D dimensional space. The positions of 
the particles on the line, Sj are determined by the equa- 
tion Ri = 0. That is, 



- 1 

g V — V i = 1, 



(20) 



1.1 



and an absolute minimum is attained when each S; sat- 
isfies this equation. 

We now show that eqs.(|2(]) in fact have a unique 
solution. We scale out g by defining a new variable, 
tj = (7 _1 / 2 Si. Define 



JV 



(21) 



Using the equilibrium condition, 
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Manipulating the summations, we get the recursion rela- 
tion, 
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Ck are clearly function of only N and are independent 
of any particular solution of eq.(p0|). Also, by definition 
C-i = N and Cq — using directly the eq. (p0|) . 

Let t'i be another solution of the eq. ( p0| ) . Because the 
Cfe are independent of the particular solution, we have, 
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where AU = U - and a kl = ^ r=0 t\{t 
plying both sides by arbitrary bk and summing over k we 
get, 
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where Pi = X]fc°=o a kibk- Since bk and therefore /3j are 
completely arbitrary we have Ati — 0. Therefore the 
solution is unique. It is well known ||, in fact, the N- 
zeros of the Hermite polynomial H^it): namely ti, satisfy 
the eq.(|20|) with ti = g~ x l 2 Si. Hence the uniqueness 
proof implies that the positions of the particles along 
a line in the ground state are determined by the zeros of 
the Hermite polynomial. 

As an aside, we also remark that since Ck are con- 
structed out of the zeros of the Hermite polynomials, we 
have also obtained an infinite set of identities for the ze- 
ros of the Hermite polynomial. 

This completes the proof that for all A > there is only 
one unique global minimum (modulo the states obtained 
by rotations). The corresponding configuration is a line 
where the positions of the particles are given as the zeros 
of Hermite polynomial in terms of the scaled variables ti. 

We have thus demonstrated that a D dimensional 
Hamiltonian with two and three body interactions has 
a unique stringlike ground state when the strength of the 
three body term exceeds a critical value determined in 
terms of the parameter A. Further we note: 

1. The result is independent of the strength of the 
two-body term, the number of particles and the 
strength of the oscillator potential. Thus even in 
the full quantum theory it is possible to choose g 
or N large enough such that the quantum fluctua- 
tions about the classical ground state are as small 
as we wish. The above statements about the classi- 
cal ground state approximately remain valid, there- 
fore, even in the quantum theory. 
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2. The string like ground state one obtains here for 
A > is identical to the classical ground state of 
the one-dimensional CSM Hamiltonian. 
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3. The model provides an example of a structural 
transition driven by the mutual interactions with- 
out any external influence. In fact this statement 
is manifestly independent of the space dimension. 
Thus we have a powerful result that in any space- 
dimension, the Hamiltonian in eq.(l) gives rise to 
string like ground state when the strength of the 
three body term exceeds the strength of the two 
body term irrespective of the number of particles, 
the strength of the confinement potential and the 
strength of the two body term. 

4. Finally, we point out that the form of the three- 
body potential as given in eq. (|l4|) , may be easily 
generalised. Note that the results of the analysis 
remain unchanged as long as V3 and its gradient 
vanish on a line and V3 is positive semi-definite. In 
particular, 

V3(rij,fjk,f k i) = (fij ■ r jk )(r jk -f k i) + cyclic 

= r 2 -rf, sin 2 (% fe ) (27) 

satisfies these properties. Hence it follows that for 
any V3 of the form, 

V3= /( r ij' r jfc' r H) w 3(^,^- fe ,f fcl ) (28) 

where /(ry, rik, fjk) is a symmetric function of its 
arguments and is positive, also vanishes for a line 
together with its gradient. The V3 considered in 
eq.(|l4|) corresponds to 

f(nj,r jk ,rki) = (^VL)" 1 

Thus, for any arbitrary / modulo these properties, 
the mere presence of a three-body interaction is 
sufficient to give a string like ground state. The 
crucial point here is that the geometric structure 
of the ground state depends on the sin 2 (%&) which 
favours string configurations and not on the precise 
form of /. The three-body interaction here differs 
from the one proposed by Weber and Stillinger || , 
in that 29ijk appears in place of 9ijk- While the 
final ground state configuration may be different 
in these two cases, the presence of the three body 
term is crucial in defining the actual ground state 
configuration. 
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